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Abstract. We consider solutions to the linear wave equation on non-compact 
Riemannian manifolds without boundary when the geodesic flow admits a 
filamentary hyperbolic trapped set. We obtain a polynomial rate of local 
energy decay with exponent depending only on the dimension. 



1. Introduction 

In this paper we consider solutions to the linear wave equation on the non- 
compact Riemannian manifolds with trapping studied by Nonncnmacher-Zworski 
jNoZw] . Let (X, g) be a Riemannian manifold of odd dimension n > 3 without 
boundary, with (non-negative) Laplace-Beltrami operator —A acting on functions. 
The Laplace-Beltrami operator is an unbounded, essentially sclf-adjoint operator 
on L 2 (X) with domain H 2 (X). We assume (X,g) is asymptotically Euclidean in 
the sense of |NoZw| (3.7)-(3.9)]. That is there exists Rq > sufficiently large that, 
on each infinite branch of M \ B(Q,Rq), the semiclassical Laplacian — h 2 A takes 
the form 



-/i 2 A| M \ B(0ii?o) = 52 a a (x,h)(hD x ) c 

\a\<2 



with a a (x,h) independent of h for \a\ = 2, 



52 a a {x,h){hD x ) a > C _1 |f| 2 , 0<C<oo, and 

a a (x, h)(hD x ) a — ► |£| 2 , as \x\ —* oo uniformly in h. 



\a\=2 

5] a 

M<2 



In order to quote the results of |NoZwj we also need the following analyticity as- 
sumption: 3#o £ [Oi 71 ") such that the a a (x,h) are extend holomorphically to 

{rw : uj e C™, dist(w,§") < e, r G C, |r| > R , argr G [-e,0 o + e)}. 

As in |NoZw| . the analyticity assumption immediately implies 



dl [ 52 a a (x,h)C~\e\ =o(\x\-^)(0 

Jo|<2 



2 , Ixl — > OO. 



We assume also that the classical resolvent (—A — A 2 ) 1 has a holomorphic 

2 r2 
comp loc ' 



continuation to a neighbourhood of A G K as a bounded operator L 2 omp — > L 2 



2 



HANS CHRISTIANSON 



We consider solutions u to the following wave equation on X x Rf. 

( (-L> 2 - A)u(x,t) = 0, (i,f)elx[0,oo) 
(l.l) < u{x,0) = u eH 1 {X)r\C™(X), 

{ D t u(x, 0) = u 1 e l 2 (x) n C£°(Jf), 

For u satisfying Ql.ip and \ € C^°(X), we define the local energy , E x (t), to be 

£ xW = ^ (\\xdtu\\ 2 L 2 (x) + \\xu\\hi(x)) ■ 

Local energy for solutions to the wave equation has been well studied in various set- 
tings. Morawetz |Mor j . Morawetz-Phillips [MoPhj ■ and Morawetz-Ralston-Strauss 
|MRS| study the wave equation in non-trapping exterior domains in R™, showing 
the local energy decays exponentially in odd dimensions n > 3, and polynomially in 
even dimensions. This has been generalized to cases with non-trapping potentials 
[Vaij and compact non-trapping perturbations of Euclidean space jVodj . In the case 
of elliptic trapped rays, it is known that (see |Ral| ) exponential decay of the local 
energy is generally not possible. Ikawa [Ikal, Ika2] shows in dimension 3 there is 
exponential local energy decay with a loss in derivatives in the presence of trapped 
rays between convex obstacles, provided the obstacles are sufficiently small and far 
apart. In the case X is Euclidean outside a compact set, dX ^ 0, and with no 
assumptions on trapping, Burq shows in |Burlj that E x (t) decays at least logarith- 
mically with some loss in derivatives. The author shows in |Chr3j that if there is 
one hyperbolic trapped orbit with no other trapping, then the local energy decays 
exponentially with a loss in derivative (including the case dX = 0). 

The main result of this paper is that if there is a hyperbolic trapped set which 
is sufficiently "thin" . then the local energy decays at least polynomially, with an 
exponent depending on the dimension n. 

Theorem 1. Suppose (X,g) satisfies the assumptions of the introduction, dimX = 
n > 3 is odd, and (X,g) admits a compact hyperbolic fractal trapped set, Ke, in the 
energy level E > with topological pressure Pg(l/2) < 0. Assume there is no other 
trapping and (—A — A 2 ) -1 admits a holomorphic continuation to a strip around 
RcC. Then for each e > and s > 0, there is a constant C > 0, depending on e, 
s, and the support of uq and U\, such that 



(1.2) E x (t) < C 



Aog(2 + i) \ 3 "+ B /„ „2 2 

y-^-j. — -J (JM&i+apo + 



Remark 1.1. It is expected that Theorem[T]is not optimal, and in fact an exponen- 
tial or sub-exponential estimate holds. Similar to in [Chr3j . we expect applications 
to the nonlinear wave equation, although there arc certain technical difficulties to 
overcome. 

The proof of Theorem [T] is a consequence of an adaptation of |Burl| Theoreme 
1] to this setting and the following resolvent estimates. 

Theorem 2. Suppose {X,g) satisfies the assumptions of Theorem]^ Then for any 
X € C^°(X) and any e > there is a constant C = C Xj£ > such that 

^-i log(l + (A)) 



\x(-A-\ z )- L x\\mx)^mx)<c- 



(A) 
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for 

AG |A. |ImA| < | C /| ReA |-3n/2-e ) | ReA |> C ]' 

Remark 1.2. The proof of Theorem [T] depends more on the neighbourhood in 
which the resolvent estimates hold than on the estimates themselves. Given a 
complex ncghbourhood of the real axis, any polynomial cutoff resolvent estimate 
will give the same local energy decay rate. Theorem [5] represents a gain over the 
estimates in |NoZwl Theorem 5] in the sense that the estimate holds in a complex 
neighbourhood of R, rather than just on ML 

Acknowledgments. This research was partially conducted during the period 
the author was employed by the Clay Mathematics Institute as a Liftoff Fellow. 



2. Proof of Theorem [2] 

To prove Theorem [2j we use the results of Nonnenmacher-Zworski [NoZw| to 
prove a high energy estimate for the resolvent with complex absorbing potential, 
then use the holomorphic continuation to bound the cutoff resolvent by a constant 
for low energies. If we consider the problem 

(2.1) (-A - \ 2 )u = /, 

and restrict our attention to values |A| > C for some constant C > 0, we can 
transform equation (|2.1[) into a semiclassical problem for fixed energy by setting 

A = y/i/h 

for z ~ 1 and < h < Hq. Then (|2 . If) becomes 

{P-z)u = h 2 f, 

where 

P = -h 2 A 

is the semiclassical Laplacian. 

The following Proposition is the high energy resolvent estimate from [NoZw 
with the improvement that the estimate holds in a larger neighbourhood of R C C. 

Proposition 2.1. Suppose W <G C°°(A"; [0, 1]), W > satisfies 

suppVK c X\ B(0,i?i), W= 1 onX\B(0,R 2 ), 

for i?2 > Ri sufficiently large, and 

( h N \ 

\\{P - iW - z)~ l \\ L ^ L 2 <C N (l + bg(l/fc) + — J , 

for z G [E — 5, E + 5] + i(—ch, ch). Then for each e > and each x € C^'(X), there 
is a constant C — C CjX > such that 

for z G [E - Cl h, E + Cl h] + i(-c 2 /i 3 "/ 2+1+e , c 2 /i 3n/2+1+e ). 
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We first improve |NoZw( Lemma 9.2] in order to get cutoff resolvent estimates 
with the absorbing potential in a polynomial neighbourhood of the real axis. The 
proof of the following lemma is an adaptation of the "three-lines" theorem from 
complex analysis and borrows techniques from |Chrll IBuZwl INoZwj and the refer- 
ences cited therein. 

Lemma 2.2. Suppose F(z) is holomorphic on 

fi = [-l,l] + i(-c_,c+), 

and satisfies 

log\F{z)\ < M, zed, 

\F(z)\ < a+— — , z G fin{Imz > 0}. 
Im z 

Then if j < eM~ 3 / 2 for e > sufficiently small, there exists a constant C = C e > 
such that 

\F{z)\ <Ca, z e [-l/2,l/2] + i(-Af- 3 / 2 ,M- 3 / 2 ). 
Proof. Choose i>(x) € C~([-l, 1]), ?/> = 1 on [-1/2, 1/2], and set 

,( 2 ) = r ./ i | e -(.-^^( x)(fc , 

where < /3 < 1 and c > will be chosen later. The function tp(z) enjoys the 
following properties: 

(a) tp(z) is holomorphic in f2, 

(b) \tp(z)\ <Confin{|Imz| < /3 1 / 2 }, 

(c) \(fi(z)\ > C^ 1 on {| Rez\ < 1/2} n {| Imz| < /?} if c > is chosen appropri- 
ately, 

(d) \ip{z)\ < Ce~ c IP for z G {±1} + z(-/? 1/2 , /3 1/2 ). 
Now for a G M to be determined, set 

g(z) = e ia *<p(z)F(z). 

For 5± > to be determined, let 

Q' :=fin{-L < Imz < 5+}. 

We have the following bounds for g(z) on the boundary of fl': 

( -C//3 + M -almz, Rez = ±1, if | Imz| < /3 1 / 2 , 

log|#(z)|<< C + M + aS-, Imz = -S->-(3 1 / 2 , 

{ C + log(a + 7 /c5 + ) -a6+, Im z = 6+ < /3 1 / 2 . 

We want to choose a, /?, and S± to optimize these inequalities. Choosing a = 
-2M/6- yields 

log|5(»| < C — M for Imz = 
and choosing 5+ = \2/a\ yields 

log \g(z)\ <C + log(a + j/5+) + 2, for Imz = 5+. 
Finally, chooing (3 = C'/M for an appropriate C > yields 

log|g(z)| < -C~ l M for Rez = ±1, |Imz| < max{<5+, <5_}, 
and taking 5- = C"M~^ 2 , 5+ = C"M~ 3 / 2 gives 

log \g(z)\ < C" + log(a + 7 /5+) on dQ'. 
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In order to conclude the stated inequality on F(z), we need to invert e laz ip(z), 
which, from the definition of a and the properties of tp stated above, is possible for 

z G [-1/2,1/2] + i(-M- 3/2 ,M~ 3/2 ). 

Then for z in this range and 7 satisfying 7 < eA/~ 3//2 , 

\F(z)\ <Ca(l + e) < C'a, 

as claimed. □ 

Now to prove Proposition l2.il as in NoZwj, we apply Lemma \2. 2 1 to 

F(Q = (h(P-iW-hC)- 1 f,g) L27 

for f,g € L 2 . For M we use the well-known estimate 

||(P- iW - z)- l \\ L *_+ L i < C 6 e ch ~ n ~\ Imz > -h/C, 

and take M = C e h~ n ~ e . For the other parameters, we take 

j = h N , a = co + log(l//i). 

Rescaling, we conclude 

\\{P- l W-z)^\\<C^^ 

in the stated region. Then we apply the remainder of the proof |NoZw| Theorem 
5]. □ 

3. Proof of Theorem [T] 

In this section we adapt the proof of [Burl! Theoreme 1] to the case where one 
has better resolvent estimates. Wc first present a general theorem on semigroups 
(see [Burl[ Theoreme 3] and jLebj ) . 

Let H be a Hilbert space, a meromorphic family of unbounded linear op- 

erators on H 7 holomorphic for Im£ < 0. Assume for Im£; < 0, 

lm(B(€)u,u) H > 0. 

Let Dom (B) = Dom(l — iB(—i)) denote the domain of B. Assume for Im£ < 0, 
£ — B(£) is bijective and bounded with respect to the natural norm on Dom(i?), 

IMI D om(B) = \nH + \\B(-i)uf H , 

and 

iKe-BCor'ik-^ciimfl- 1 . 

Assume that B(£) G 5 1 (]R 2 ; £(Dom (B),H)). That is, B(£) is a symbol with respect 
to £ and assume that, as operators on Dom (B). 

lU !):„■< = e is tB(£ + D s ) 

= e is ^B(0, 

since members of Dom (B) do not depend on s. We assume B satisfies the identity 

B(D t )tl>(t)U(t) = i>(t)B(il>'/iip + D t )U(t) 
= m(B(D t ) + A(t))U(t), 

for ip(t) e C°°(M), and U G C°°(R t ; Dom (B)). Here, A(t) is a linear operator, 
bounded on H and has compact support contained in supp^'- 
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By the Hille-Yosida Theorem, for every k £ N and s > 0, we can construct the 
operators 

e isB{D s ) 



0= (eC: |Im£| < 



(1 -iB{-i)) k ' 

where e lsB ^ D "' satisfies the evolution equation 

(D s - B(D s ))e isB ^ = 0, 
e-^)| s=0 = id. 

Now suppose Xji 3 = 1, 2 are bounded operators i? — > iJ, and xi(£ — ^(£)) _1 X2 
continues holomorphically to the region 

G, |Rc£|<G 
P(|Re£|), |Re£|>G, 

where P(|Re£|) > and is monotone decreasing (or constant) as |Re£| — > oo. 
Assume 

(3.1) ||xi(£ - B(0)- 1 X2\\h^h < G(\ Re^l) 

for £ e fi, where G(| Re£|) = Q(\ Re£\ N ) for some N > 0. We further assume that 
the propagator e lsB ^ DB ^ "acts finitely locally," in the sense that for s € [0, 1], 

is also a bounded operator on H, and Xi(£ — P(£)) _1 X2 continues holomorphically 
to fl and satisfies the estimate (|3.1[) with G replaced by GG for a constant G > 0. 



Theorem 3. Suppose satisfies all the assumptions above, and let k G N, 

fc > A + 1. T/ien /or any P(i) > ; monotone increasing, satisfying 

(3.2) P(t) fe+1 < exp(fP(F(f))), 

we have 

e ttB{D t ) 



(3.3) 



Xl (l -iB(-i)) kX2 



< CF(ty 



As in [Burl] . Theorem [T] follows from Theorem [3] by setting 
B - ^ -iA 

the Hilbert space if = H\X) x L 2 (A), and X j e C C °°(A) for j = 1,2. The 
commutator [x2,P] is compactly supported and bounded on H, so if X2 G C£°(X) 
is supported on a slightly larger set than X2, we have 

llxie JtB X2||Do m( B^ = \\xie UB X2(l-iB)- k \\ H ^ H 

< C\\ X ie itB (l-iB)- k x 2 \\ M . 
Taking jfe = 2, P(t) = t-3n/2-e/2 ) and 

yields (|1.2[) for s > k. We observe the spaces i? 1+s x fP are complex interpolation 
spaces, hence interpolating with the trivial estimate 

E x (t)< \\uo\\h + hi\\l^ 
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yields (JOJ) for s > 0. 

□ 

Remark 3.1. Evidently, if we have polynomial resolvent bounds in a fixed strip 
around the real axis, we have exponential local energy decay for the wave equation 
with a loss in derivatives. Further, if H = L 2 (X) for X a compact manifold, 
this theorem may be applied to the damped wave operator with \i = X2 — 1 to 
conclude there is exponential energy decay with loss in derivatives for solutions to 
the damped wave equation if there is a polynomial bound on the inverse of the 
damped wave operator in a strip. This corrects a mistake in the proof of [Chrl| 
Theorem 5]. 

We first need a lemma. 

Lemma 3.2. For k > N + I, the propagator satisfies the following identity on H : 

p itB(D t ) 1 f 

(\-iB(-i)) k 2mJ Im ^_i 

Proof. We write Ik for the right hand side and observe both the left hand side and 
Ik satisfy the evolution equation 

(A - B(D t ))w = 0. 

To calculate ifc(0), we deform the contour to see 

4(0) = f / - / ] (l - i0~ k ^ - Btf))- 1 *. 

Letting C — > oo, the first integral vanishes. Thus we need to calculate the second 
integral. For k = 1, this is the residue formula, while for k > 1 the formula follows 
by induction and the continuity of B(£) as e — > 0. 

Thus the left hand side and Ik have the same initial conditions, and the lemma 
is proved. □ 

Proof of Theorem [3J Now, as in |Burl| , we introduce a cutoff in time to make the 
equation inhomogeneous, then analyze the integral separately for low and high 
frequencies in £. In order to maintain smoothness, we convolve with a Gaussian. 
For an initial condition uq G H, let V(t) = e ltB ^ Dt ' , X2UQ^ and consider U(t) = 
i/)(t)(l - iB(-i))- k V(t) for ip(t) G C°°(R) satisfying ip = for t < 1/3, ip = 1 for 
i > 2/3, and V' > 0. We observe by the sub-unitarity of e ltBi - D ^ for t > 0, 

\\U(t)\\ <C\\V(t)\\ <C'\\u \\, 

where for the remainder of the proof, | ■ || = || ■ \\h unless otherwise specified. 
The family U{t) satisfies 

(A - B{D t ))U = A(t)(l - iB{-i)y k V{t), 

where A is a bounded operator on H with support contained in [1/3,2/3]. As 
U(0) = 0, Duhamel's formula yields 

U(t)= f e l{t - s)B{Dt) A(s){l-iB{-i))- k V(s)ds, 
Jo 
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and by Lemma 13. 2 



U(t)= / e« t -'KA(s)(l-it)- k (Z-B(0)- 1 V{s)dZda. 

Js=0 J Im£=-l/2 

For a function F(i) > 0, monotone increasing in t to be selected later, we will cut off 
frequencies in |£| above and below F{t) 2 . We convolve with a Gaussian to smooth 
this out: 

U(t) = f [ /e'C'-'KAMl -A®)" 1 

Js=0 J Imf=-l/2 J\ 

{c /TT)h- C0 ^-^ F ^ 2 V(s)dXd^ds 

= f [ ( [ + [ J (-)dA^ S 

JO ./Im£=-l/2 W|A|<-F(t) J\\\>F(t)J 
= : h+h- 

Analysis of I-^: From the resolvent and propagator continuation properties, the 
integrand in I\ is holomorphic in { Im£ < 0} U il. Observe if | Re£| 3> F(t) 2 , then 
the integrand is rapidly decaying, hence we can deform the contour in £ to 

C, | Rc^l < C 

P(|Re£|), |Re£|>C. 

We further break I± into integrals where Re£ is larger than or smaller than 
F 2 (t): 

h = f I [ + f ) I (-)dXd^ds 

JO \Jrn{\ Re£|<AF(t) 2 } JTn{\ Re£|>AF(i) 2 } J J\\\<F(t) 
= : Jl + J 2 - 

For Ji, if t > 2, since P(| Re£|) is monotone decreasing, we have 

Im£ > P(AF(t) 2 ), 
and on the support of A, we have t — s > t — 1. Hence 



r = U e C : Im£ = 



C 



/ / 

Jrn{|ReC|<A_F(t) 2 } J|A|<F(t) 



r (*-D^(*) J )^)-* G (|He^|) 



-c (A-£/F(t))< 



dAde||uo| 



For J2, we observe that for A large enough and | Re£| > AF(t) 2 , 
Re (A- £/F(t)) 2 > C~ l (X 2 + (Re® 2 /F(t) 2 ). 



Hence, 



\X1J2 



C 



C 



rn{| RoC|>A_F(t) 2 } J\x\<F(t) 



(0- k G(\Ret\) 



-c (\^/F(t)Y 



dXd£\\v.Q\ 



F(t)e 



'\n\>F(t) 
< CF(t)e- c ^\\u \\. 



dr)\\u \\ 
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Analysis of li: Set 



J(r) = I _ A(s)e^-^^-iO- k ^-B(0)- 1 



■(c /7r)h- Co ^-^ F ^ )2 V(s)dXd^ds, 
which for t > 1 is equal to U (r) . Observe 

(D T - B(D T ))J(r) = f f A(s)e^- S ^(l - 

Js=0J j^|> F(t) 

■(c /TT)h~ C0 ^ /F ^ 2 V{ S )dXd^ds 

= ■■ K(t). 

Hence 

J(t) = e ltB{D '\7(0) + /' e^- s)B ^K(s)d S . 



Again, by the subunitarity of the propagator, we need to estimate ||J(0)|| and 
Jp ||-K"(s)||ds. For s <G since k > N + 1, we can deform the ^-contour in the 

definition of K to Im£ = F(t). Then for this range of s, 

\\K(s)\\<C [ e-(- 2 /3)F( t ) (^-^Huoll, 

and hence 



J \\K(s)\\ds <CF(ty 



l e -F(t)/3_ 



For J(0), we first consider A > F(t). Since k > N + 1, we can deform the 
^-contour to 

r' = r_ u r+ 

where 

T_ = {Rc£ < F(t) 2 /A, Im£ = -1/2} 

U{ Rc£ = F(t) 2 /A, < Im < -1/2} 

and 

T + = {Ret >F(t) 2 /A, lm£ = -F(t)}. 

If £ e r_, we have 

Rc(A-e/^(t)) 2 > X 2 /C, 

so 

(0" fe G(|ReC|) • e- C0 ^/ F ^ 2 V(s)d\dt < Ce~ F ^ . 



l e -istl = e -F(t)/3 



/?er_ J\>F(t) 
For £ € r + , we have 



so the contribution to || J(0)|| coming from A > F(t) is bounded by 

C(e^ 2 + e - F «/ 3 ). 
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The contribution to ||J(0)|| coming from A < —F(t) is handled similarly to obtain 
the same bound. 

Wc have yet to estimate L \\K(s)\\ds. For this we use Plancherel's formula to 
write 



\K(s)\\ds 



< 



K{s)\\ 2 ds 

(i-i€)-*2v(o 



(3.4) = / (l-%)- k AV($) e -«=o(A-€/F(t)) dX ^ 

J-oo J\\\>F(t) 

If we estimate this integral by again considering regions where |£| < F(t) 2 /A and 
\£,\ > F(t) 2 /A respectively, we see (|3.4j) is majorized by 



C(F(t) 



-2k 



-F(ty/c 



C{F{ty 2k +e~ F ^ 2 / c ) / \\AV{s)fds 



< c{F(ty 2k + e- F ^ 2 ' c )\\u4 2 . 

Combining all of the above estimates, we have 

F{t)~ 2k 



\\U{t)\\ < Cmax 



-F(t)/3 



-F(ty/c 



\ u o\ 



F(t)- 1 e- F «/ 3 , 
F{tfe- tP ^tf) + F(t)e F ® 

Relabelling F(t) 2 as F(t) throughout and applying the condition (|3.2p . we recover 



□ 
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